Properties of light resonances from unitarized Chiral perturbation 
theory: Nc behavior and quark mass dependence 
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We review the unitarization of Chiral Perturbation Theory with dispersion relations and 
how it describes meson-meson scattering data, generating light resonances whose mass, width 
and nature can be related to QCD parameters like quark masses and the number of colors. 

§1. Introduction 

o' 

CN ■ Light hadron spectroscopy lies beyond the apphcabihty of perturbative QCD. How- 

ever, there is an effective field theory, known as Chiral Perturbation Theor}© (ChPT), 
which provides a description of the dynamics of the lightest mesons. Despite it is lim- 
f-^ • ited to low energies and masses, here we review how, when combined with dispersion 

^SJ I relations, it leads to a successful description of meson dynamics, generating resonant 

states without a priori assumptions on their existence or nature. This "unitarized 
,i^ . ChPT" is a useful tool to identify the spectroscopic nature of resonances through 

Mh| their dependence on the QCD number of colors N^, but also to relate lattice results 

^1^' to physical resonances by studying their quark mass, niq, dependence. 

ChPT is built out of the Goldstone Bosons of the QCD spontaneous chiral symmetry 
breaking, namely, pions, kaons and etas, as a low energy expansion of a Lagrangian 
respecting all QCD symmetries. It is organized in powers oi p'^/A'^, where p stands 
either for derivatives, momenta or meson masses, and A = 47r/^, where /tt denotes 
the pion decay constant. ChPT is renormalized order by order by absorbing loop di- 
22 ' vergences in the renormalization of parameters of higher order counterterms, known 

^f-N . as low energy constants (LECs) that carry no energy or mass dependence and de- 

pend on a regularization scale /x. As always after renormalization, the full amplitude 
t^ ' is independent of this scale. Their values depend on the QCD dynamics, and are 

^D . determined from experiment. Up to the desired order, the ChPT expansion provides 

a systematic and model independent description of how meson observables depend 
on QCD parameters like the light quark masses m = {niu + rnd)/'2 and tti^, or the 
^ . leading \/Nc behavior .^^ 

Cd \ §2. Dispersion relations and unitarization 

Elastic resonances appear as poles on the second Riemann sheet of the meson- 
meson scattering partial waves t/j of definite isospin / and angular momentum J . 
At physical values of s, elastic unitarity implies 

Imt7j(s) = a(s)|t/j(s)|2 ^ tij = --—^ -, with a{s) = 2p/^s, (2-1) 

Re tjj — la 
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where s is the Mandelstam variable and p is the center of mass momentum. However, 
ChPT amplitudes, being an expansion t/j ~ t\j + t\j + • • • , with i*^^*^) = 0{p'^^), 
can only satisfy Eq. ()2-ip perturbatively 



Im45(.) = 0, lu.tfj{s) = a{s)\tu{si'^\ •••^ Im4^j(.)/4f (.) = ct(.), (2-2) 

and cannot generate poles. However, the resonance region can be reached combining 
ChPT with dispersion theory either for the amplitude^^ or for the inverse amplitude 
through the Inverse Amplitude Method (IAM).'^t^ We will concentrate on the one- 
channel IAM,™'I3 since it uses ChPT only up to a given order inside a dispersion 
relation, without additional input or further model dependent assumptions. Other 
unitarization techniques will be commented below. 

2.1. The one-loop ChPT Inverse Amplitude Method 

For a partial wave tij{s), we can write a dispersion relation (that we subtract three 
times, since we will also use it below for tjj, that grows with s^) 

tuis) = Co + C,s + C^s' + - r ]!^l''^''^'^'\ + LC{tu). (2-3) 

Note we have explicitly written the integral over the physical cut, extending from 
threshold, sth^ to infinity, but we have abbreviated by LC the equivalent expression 
for the left cut (from to — oo). We could do similarly with other cuts, if present, as 
for TiK — >• ttK. Note that from Eq. ()2Tp the imaginary part of the inverse amplitude 
is exactly known in the elastic regime. We can then write a dispersion relation like 
that in ()2-3p but now for the auxiliary function G = {tjj)'^ /tij, i-e., 

G{s) = Go + Gis + G2S^ + - r ^^^^''^"^''^ + LG{G) + PG, 

where now PC stands for possible pole contributions in G coming from zeros in tjj. 

(2) 

It is now straightforward to expand the subtraction constants and use that Im tjJ = 
and Imtjj = (t|4jP, so that ImG = — Imt^-j. In addition, up to the given order, 
LG{G) ~ —LG{tjj), whereas PC is of higher order and can be neglected. Then 

^^.„,... ._ ... U^.2 S^ r lmtf]{s')ds' ,^,,(4), ^,(2) ,(4) 

'Sth 



— _„„ + „,.- 6„ - M - M^ --/ jsij^ifji-j - Lcw;) =. ty - ti, . 



since the Oj, bi terms, coming from the Gi expansion, are the subtraction terms of a 
dispersion relation for tjj — tjj . Thus we arrive at the so-called lAM: 

t/.-4f/(45-«), (2-4) 

that provides an elastic amplitude satisfying unitarity and has the correct ChPT 
expansion up to the order we have used. The PG contribution has been calculated 
explicitly™' and is not just formally suppressed, but numerically negligible except 
near the Adler zeros, away from the physical region. It is straightforward to extend 
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the lAM to other elastic channels or higher orders.'^ Naively, the I AM looks like 
replacing Ret^j by its 0{p^) ChPT expansion in (j2-ip . but ()2Tp is only valid in the 
real axis, whereas our derivation allows us to consider the amplitude in the complex 
plane and look for poles associated to resonances. Let us remark that, since ChPT 
is used only at low energies in the dispersion relation, the lAM formula is justified 
only up to energies where inelasticities become important, even though ChPT does 
not converge at those energies. Only when the energy is close to the Adler zero one 
should use a slightly modified version of the lAM.® Re-expanding the lAM, ChPT 
is recovered up to the order it was used as input, as well as partial contributions to 
higher order, but not the complete series — see Refl?]) for a discussion of this issue. 
In Fig d] we present some resultsS of an updated fit of the lAM tttt and ttK scatter- 
ing amplitudes to data, simultaneously fitting the available lattice results on ratios 
of meson masses and decay constants and some scattering lengths. It is important 
to remark that the resulting LECs are in fairly good agreement with standard deter- 
minations: no fine tuning is required. The /o(600), /9(770), k(800) and i^*(892) are 
not introduced by hand but generated as poles in the second Riemann sheet of their 
corresponding partial waves. The fact that we do not need to model the integrands 
and the only input parameters are those of ChPT is relevant since we then know how 
to relate our amplitudes to QCD parameters like Nc or niq. 

2.2. Other unitarization techniques within the coupled channel formalism 

Naively one can arrive at ()2-4p in a matrix form, ensuring coupled channel unitarity, 
just by expanding the real part of the inverse T matrix. Unfortunately, there is still 
no dispersive derivation including a left cut for the coupled channel case. Being much 
more complicated, different approximations to ReT~^ have been used: 

• The fully renormalized one-loop ChPT calculation of ReT"^ provides the correct 
ChPT expansion, with left cuts approximated to O(p^).'^''^ Indeed, using LECs 
consistent with standard ChPT determinations, one can describelSI' below 1.2 GeV 
all two-body scattering channels made of pions, kaons or etas. Simultaneously, this 
approach 1^ generates poles associated to the p(770) and K*{892) vector mesons, 
together with the /o(980), ao(980), /o(600) and k (or Kq{800)) scalar resonances. 

• Originally, ^'^■' the coupled channel lAM was used neglecting crossed loops and 
tadpoles. This is considerably simpler, and despite the left cut is absent, since its 
numerical influence is relatively small, meson-meson data are described with reason- 
able LECs while generating all poles enumerated above. Note that this approxima- 
tion keeps the s-channel loops but also the tree level up to O(p^), which encodes the 
effect of heavier resonances, like the p. Thus, contrary to some common belief, the 
lAM incorporates the low energy effects of t-channel p exchange. 

• Finally, if only scalar meson-meson scattering is of interest, it is possible to use 
just one cutoff (or another regulator) that numerically mimics the combination of 
LECs appearing in scalar channels. This " 'chiral unitary approach" is very popular, 
even beyond the meson-meson framework, due to its great simplicity but remarkable 
success^ and also for its straightforward relation to the Bethe-Salpeter formalism^^ 
that provides physical insight on unitarization. With this method it was showipSJ 
that, assuming no niq dependence of the cutoff, all light scalar resonances degenerate 
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Fig. 1. Updated lAM fit*' (continuous line). We also show non-unitarized ChPT results with the 
LECs from the Ku two-loop analysis^ (dot-dashed line). Left: lAM versus data on tttt and -kK 
scattering. Right: fit results compared to lattice calculations^^- on ratios of meson masses and 
decay constants and some scattering lengths. We fit up to 7n^ = 440 MeV, but even beyond 
(grey areas) lattice results are not described badly. Experimental references are detailed in llOp . 

into an octet and a singlet in the SU(3) limit. Axial- vector mesons have also been 
generated by using a chiral Lagrangian for the pseudoscalar-vector interaction.'^ 

§3. The nature of resonances from their leading 1/Nc behavior 

The QCD 1/Nc expansion^ valid in the whole energy region, provides a rigorous 
definition of qq bound states: their masses and widths behave as 0(1) and 0{1/Nc), 
respectively. The QCD leading 1/Nc behavior of Jt^ and the LECs is well known, 
and ChPT amplitudes have no cutoffs or subtraction constants where spurious Nc 
dependences could hide. Hence, by scaling with Nc the ChPT parameters in the lAM, 
the Nc dependence of the mass and width of the resonances has been determined to 
one and two loops.'^''^ These are defined from the pole position as ^Spde = M—iF. 
However, a priori, one should be careful not to take Nc too large, because the Nc ^- oo 
limit is a weakly interacting limit. As shown above, the lAM relies on the fact 
that the exact elastic RC contribution dominates the dispersion relation. Since the 
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Fig. 2. Left: p (top) and a (bottom) pole trajectories for different values of /i, note that for p, = 1.2 
GeV the p pole goes away the real axis. Center: Nc behavior of the p (top) and K* (bottom) 
mass and width. Right: Nc behavior of the a mass and width. 

lAM describes data and the resonances within, say, 10 to 20% errors, this means 
that at A'^c = 3 the other contributions are not approximated badly. But meson 
loops, responsible for the RC, scale as 2>/Nc whereas the inaccuracies due to the 
approximations scale partly as 0(1). Thus, we can estimate that those 10 to 20% 
errors at A'^c = 3 become 100% errors at, say Nc ~ 30 or Nc ~ 15, respectively. Hence 
we never show results^''^ beyond Nc = 30. Even beyond Nc ~ 15 they should be 
interpreted with care. 

Thus, Figi2] shows the behavior of the p, K* and a masses and widths found va.)^ 
The p and K* neatly follow the expected behavior for a qq state: M ~ 1, I^ ~ ^/Nc- 
The bands cover the uncertainty ^ ~ 0.5 — 1 GeV where to apply the 1/Nc scaling. 
Note also that outside this p range the p meson starts deviating from a qq behavior. 
Something similar occurs to the A'*(892). Hence, we cannot apply the Nc scaling at 
an arbitrary p value, if the well established p and K* qq nature is to be reproduced. 

In contrast, the a shows a different behavior from that of a pure qq: near Nc=3 
both its mass and width grow with Nc, i.e. its pole moves away from the real axis. 
Of course, far from Nc = 3, and for some choices of LECs and p, the a pole might 
turn back to the real axis,'^"'^ as seen in Fig. 2 (top-right). But, as commented 
above, the JAM is less reliable for large Ac, and at most this behavior only suggests 
that there might be a subdominant qq component. ^^' In addition, we have to ensure 
that the LECs fit data and reproduce the vector qq behavior. 

Since loops are important in determining the scalar pole position, but are 1/Nc 
suppressed compared to tree level terms with LECs, we checked the O(p^) results 
with an 0{p^) lAM calculation in SU{2)W^ We defined a x^-like function to measure 
how close a resonance is from a qq Nc behavior. First, we used it at 0{p'^) to show 
that it is not possible for the a to behave predominantly as a qq while describing 
simultaneously the data and the p qq behavior, thus confirming the robustness of 
the conclusions for Nc close to 3. Next, we obtained a 0{p^) data fit where the 
p qq behavior was imposed (see Fig.3, left and center). Note that both M„ and 
Fa grow with A^r; near A^^ = 3, confirming the O(p^) result of a non qq dominant 
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Fig. 3. Left and center: Nc behavior of the p and a pole at 0{p^) with the "p as qq fit". Right: 
Sigma behavior with Nc at 0{p^) with the "p and a as qq fit". 

component. However, for Nc between 8 and 15, where we still trust the lAM, Mo- 
becomes constant and F^j starts decreasing. This may hint to a subdominant qq 
component, arising as loops become suppressed when Nc grows. Finally, by forcing 
the a to behave as a qq, we found that in the best case (Fig.3, right) this subdominant 
qq component could become dominant around A'^c > 6 — 8, at best, but always with 
an Nc ^ oo mass above ~ 1 GeV instead of its physical ~ 450 MeV value. 

Let us emphasize agairP^ what can and what cannot be concluded from our results 
and clarify some frequent questions and doubts: 

• Most likely, scalars are a mixture of different states, but the dominant component 
of the a and k in m,eson-m,eson scattering does not behave as a qq. If the qq was 
dominant, they would behave as the p or the K* in Fig. 2. However, a smaller 
fraction of qq cannot be excluded and is somewhat favored in our 0{p^) analysis.l^ 

• Two meson and some tetraquark states^ have a consistent "qualitative" behavior, 
i.e., both disappear in the meson- meson scattering continuum as Nc increases. Our 
results are not able yet to establish the nature of that dominant component. To do 
so other tools^''^ might be necessary. The most we can state is that the behavior 
of two-meson states or some tetraquarks might be qualitatively consistent. 

The Nc ^ oo limit has been studied in [20 |) . l2ip . Apart from its mathematical 
interest, it could have some physical relevance if the data and the large Nc uncertainty 
on the choice of scale were more accurate. Nevertheless: 

• A priori the lAM is not reliable in the Nc ^- oo limit, since that is a weakly 
interacting theory, where exact unitarity becomes less relevant in confront of other 
approximations made in the lAM derivation. It has been showrp^l) that it might 
work well in that limit in the vector channel of QCD but not in the scalar channel. 

• Another reason to keep Nc not too far from 3 is that we have not included 
the 7/'(980), whose mass is related to the Ua{^) anomaly and scales as y^S/Ac- 
Nevertheless, if in our calculations we keep Ac < 30, its mass would be > 310 MeV 
and thus pions are still the only relevant degrees of freedom in the a region. 

• Contrary to the leading 1/Nc behavior in the vicinity of Nc = 3, the Nc ^ oo limit 
does not give information on the "dominant component" of light scalars. The reason 
was commented above: in contrast to qq states, that become bound, two-meson and 
some tetraquark states dissolve in the continuum as Nc — )• oo. Thus, even if we 
started with an infinitesimal qq component in a resonance, for a sufficiently large Nc 
it may become dominant, and beyond that A^c the associated pole would behave as 
a qq state. Also, since the mixings of different components could change with A'c, a 
too large Nc could alter significantly the original mixings. 
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Actually, this is what happens for the one- loop lAM a resonance for Nc ^ oo, but 
it does not necessarily mean that the "correct interpretation [...] is that the a pole 
is a conventional qq meson environed by heavy pion clouds". ^^' That the a is not 
conventional is simply seen by comparing it with the "conventional" p and K* in 
Fig. 2. A large two-meson component is consistent, but so is a tetraquark. Actually, 
the A/'c — )• oo of the one-loop unitarized ChPT pole in the scalar channel limit is not 
unique^''^ given the uncertainty in the chiral parameters. Moreover, despite the 
one-loop lAM could make sense in the Nc ^ oo limit for the vector channel,'22J' in 
the scalar channel it can lead to phenomenological inconsistencies^ for some LECs, 
since poles can even move to negative squared mass values (weird), to infinity or 
to a positive mass square. Hence, robust conclusions on the dominant light scalar 
component can be obtained not too far from real life, say Nc < 15 or 30, for a /i choice 
between roughly 0.5 and 1 GeV, that simultaneously ensures the qq dependence for 
the p and K* mesons. Note, however, that under these same conditions the two- 
loop lAM still finds, not only a dominant non-qq component, but also a hint of a qq 
subdominant component, ^^' which is not conventional in the sense that it appears 
at a much higher mass than the physical a. This subdominant component at that 
higher mass seems to be needed to ensure fulfillment of local duality'^ili for Nc > 3. 
This may support the existence of a second scalar octet, a qq now, above 1 GeV^ 

Finally, using not the I AM, but the chiral unitary approach with a natural range 
for the cutoff Nc dependence, it has also been suggested—'* that a large, in some 
cases dominant, non qq behavior could exist in axial vector mesons. 

§4. Quark mass dependence of resonances 

ChPT provides a rigorous expansion of meson masses in terms of niq (at leading 
order M'^cson ~ "*-<?)• Thus, by changing the meson masses in the amplitudes, we 
see how the poles generated with the lAM depend on mq. We report here the SU(2) 
analysis^ of p and a as well the SU(3) analysis- ' of non-strange, p and a, and 
strange, k(800) and i^*(892), resonances. 

The values of m.,^ considered should fall within the ChPT range of applicability and 
allow for some elastic vrvr and ttK regime below KK or Kr] thresholds, respectively. 
Both criteria are satisfied if m^ < 440 MeV, since SU{2>) ChPT still works with 
such kaon masses, and because for m.,^ ~ 440 MeV, the kaon mass becomes ~ 600 
MeV. Of course, we expect higher order corrections, which are not considered here, 
to become more relevant as m-,^ is increased. Thus, our results become less reliable 
as m^ increases due to the 0{p^) corrections which we have neglected. 

Fig. m (left) shows the evolution of the a and p pole positions as m.,^ is increased. 
In order to see the pole movements relative to the two pion threshold, which is also 
increasing, we use units of m-,^^ so the threshold is fixed at ^/s = 2. Both poles 
move closer to threshold and they approach the real axis. The p poles reach the real 
axis at the same time that they cross threshold. One of them jumps into the first 
sheet and stays below threshold in the real axis as a bound state, while its conjugate 
partner remains on the second sheet practically at the very same position as that in 
the first. In contrast, the a poles go below threshold with a finite imaginary part 
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Fig. 4. Left: Movement of the o (dashed hnes) and p (dotted hnes) poles for increasing m.n- 
(direction indicated by the arrows) on the second sheet. The filled (open) boxes denote the pole 
positions for the a (p) at pion masses m^ = 1, 2, and 3 x m^^^ , respectively. For m^ — Zm^^^ 
three poles accumulate very near the threshold. All poles are always far enough from the Adler 
zero (circles). Right: Comparison of our results for the Mp dependence on m-,r with some recent 
lattice results.'23 The grey band covers the error coming from the LECs uncertainties. 

before they meet in the real axis, still on the second sheet, becoming virtual states. 
As rriTj increases, one pole moves toward threshold and jumps through the branch 
point to the first sheet staying in the real axis below threshold, very close to it as m.,^ 
keeps growing. The other a pole moves down in energies away from threshold and 
remains on the second sheet. These very asymmetric poles could signal a prominent 
molecular component ,123 'ESI ^t least for large pion masses. Similar movements were 
found within quark models^'^' and a finite density analysis.'^ 

Fig. S] (right) shows our results for the p mass dependence on m.^- compared with 
some recent lattice results,"^ and the PDG value for the p mass. Now the mass is 
defined as the point where the phase shift crosses 7r/2, except for those m,^ values 
where the p becomes a bound state, where it is defined again from the pole position. 
Taking into account the incompatibilities between different lattice collaborations, 
we find a qualitative good agreement with lattice results. Note also that the ttit^ 
dependence in our approach is correct only up to NLO in ChPT, and we expect higher 
order corrections to be important for large pion masses. The Mp dependence on m,r 
agrees also with estimations for the two first coefficients of its chiral expansion.l^S 

In Fig. (left) we compare the m-.^- dependence of Mp and M^- (defined from the 
pole position ^Spoie = M — ir/2), normalized to their physical values. The bands 
cover the LECs uncertainties. Both masses grow with 771,^, but M„ grows faster 
than Mp. Below m,7r ~ 2.4 m^ ^^ we only show one line because the two conjugate 
a poles have the same mass. Above 2.4 m,^ ^^, these two poles lie on the real axis 
with two different masses. The heavier pole goes towards threshold and around 
m,r — 3.3 m^ ^^ moves into the first sheet, but that is beyond our applicability limit. 

In the next panel of Fig. [5] we compare the m-j^ dependence of Fp and Fu normalized 
to their physical values: note that both widths become smaller. We compare this 
decrease with the expected phase space reduction as resonances approach the tttt 
threshold. We find that Fp follows very well this expected behavior, which implies 
that the pvrvr coupling is almost m^ independent. In contrast, F^ deviates from the 
phase space reduction expectation. This suggests a strong m.^- dependence of the a 
coupling to two pions, necessarily present for molecular states.'22J'''23li 
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Fig. 5. m-^ dependence of resonance masses and widths in units of the physical values. In the two 
left panels the dark (light) band shows the results for the a (p). The band width reflects the 
uncertainties in the SU(2) LECs. Similarly, the two right panels, calculated within SU(3),'^ 
show the behavior for the 7^* (892) (continuous) and k(800) (dashed). The (dotted) dot-dashed 
line shows the m^ dependence of the corresponding vector (scalar) width from the change of 
phase space only, assuming a constant coupling of the resonance to two mesons. 

Finally, in the last two panels of Figl5]we compare the mass and width dependence 
on rh of the ac(800) versus the i^*(892), keeping nis fixed.'S'' Note that the same 
pattern of the a — p system is repeated. Belonging to the same octet, K*{892) and p 
behave very similarly, and both their widths follow just phase space reduction. The 
a and k behaviors are only qualitatively similar, the latter being somewhat softer. 
This might be partly due to a possible significant admixture of singlet state in the 
a. The dependence of these resonances on nig has been also studied in RefEj). 

§5. Summary 

We have reviewed how the Inverse Amplitude Method (lAM}^ is derived from the 
first principles of analyticity, unitarity, and Chiral Perturbation Theory (ChPT) at 
low energies. It is able to generate, as poles in the amplitudes, the light resonances 
appearing in meson-meson elastic scattering, without any a priori assumptions. Up 
to a given order in ChPT, it yields the correct dependences on rh, rag and A'^c- 

The leading l/N^ behavior suggests that the dominant component of light scalars 
does not behave as a qq state as Nc increases not far from N^ = 3. When using the 
two loop lAM result in SU(2), below Nc ~ 15 or 30, there is a hint of a subdominant 
qq component, but arising at roughly twice the mass of the physical a. 

We have studied the pion (quark) mass dependence of the /o(600),p(770), k(800) 
and i^*(892) poleaSI''EII' and how they become bound states: softly for vectors and 
with a non-analyticity for scalars. We found that the vector-meson-meson coupling 
constant is almost m.,^ independent and a qualitative agreement with some lattice 
results for the p mass evolution with m^. These results may be relevant for studies 
of the meson spectrum^^' and form factors^ on the lattice. 
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